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A number of new properties of the spheroidal functions of arbitrary real order a > — 1 are estab- 
lished, including double orthogonality over two separate intervals simultaneously and the existence 
of a new kind of characteristic numers y««(c) that arise from it. Some computational formulas are 
derived and a few numerical results are shown. 
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1 . Introduction 

The spheroidal functions as defined and investigated by Stratum [1] ! and later by Chu and 
Stratton [2] are those solutions of the differential equation 

(l-7? 2 )l#in(c, T7)-2(a+l)7Jl/fi„(c, T?) + (6 a „-cV)^an(c, 7?)=0 (1) 

that remain finite at the singular points r\ =± 1. The condition of finiteness restricts the admissible 
values of the parameter b a n(c) to a discrete set of eigenvalues, indexed by n = 0, 1, 2, 3, . . ., that 
depend upon the values chosen for the order a and the parameter c. The differential equation (1) 
together with the condition of finiteness at T)=±l is equivalent to the integral equation [1(27)] 



Van(c)l}fan(c, T)) = j ' e ic ^(l ~ t 2 ) «l/f„ M (c, t)dt, - 1 ^ T) ^ 1, 



(2) 



with eigenvalues v an {c), valid for real a > — 1. 

For integral orders a = m = 0, 1, 2, 3, . . . the functions i// a/ ,(c, 77) are related to the spheroidal 
wave functions S m ,m+n(c, 7)) that arise from separation of the wave equation in spheroidal co- 
ordinates. For half-integral orders a=±i they are related to the periodic Mathieu functions 
Se n (c, 77) and So n + \(c, 17) that arise from separation of the wave equation in elliptic cylinder 
coordinates. The purpose of the present paper is to derive some new properties of the spheroidal 
functions of arbitrary real order a > — 1, including a new kind of characteristic numbers y a n(c) that 
arise from the property of double orthogonality shown here to be possessed by these functions, and 
to present some numerical results illustrating them. It will be assumed throughout that the param- 
eter c is real and positive, although many of the results are more general. A secondary purpose is 
to emphasize the underlying unity of the properties associated with the spheroidal functions of 
arbitrary real order a > — 1. All properties of the spheroidal wave functions and of the periodic 
Mathieu functions arise quite naturally from those of the general functions ijj a n(c, tj). An extensive 
and completely general treatment for all solutions of (1) for arbitrary complex values of a, 6, and c 
has been developed by Meixner and Schafke [3]. 
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2. Elementary Properties 

The following properties of the spheroidal functions \\f a n{c,r\) are evident directly from either 
the differential equation (1) and its boundary condition of finiteness at tj =±1 or the equivalent 
integral equation (2). 

(a) The functions i//«„(c, 77) are entire functions of 77. This follows from the fact that the param- 
eter ban was chosen in such a way that one of the two solutions of differential equation (1) remains 
finite at the singular points 17 = ±1. With no singularities except for the one at infinity that solution 
tyan(c, 77) must be entire. 

(b) As c-*0 the functions i|/«/*(c, 17) become proportional to the Gegenbauer functions 7^(77), 
and the eigenvalues b an {c) approach 

ba ti (c) -> n(n + 2a -hi). (3) 



This follows from the fact that the differential equation (1) approaches the Gegenbauer differential 
equation [4, p. 783] 

(l-T ? 2 )rf(7 ? )-2(a+l)r ? r« , (7 ? ) + Ai(^ + 2a+l)r«(7 ? )=0. (4) 

The Gegenbauer functions for a = m = 0, 1, 2, 3, . . . are related to the associated Legendre 
functions P™+n (v) DV 

1 n\VJ q_ 2 )m/2' W 

and for a = — ^ and +i they are related to the Chebyshev polynomials of the first and second 
kinds, respectively, by 

nT l ~Xiq)=J — cos [rccos^Tj], (6) 

2 sin [(n-\- 1) cos ^hj] 
77 VT^V 



y;w=v- vi-- 2 (7) 



(c) For c— >0 and 77 ^ ^ such that C77 remains finite the function *//«„(<:, 77) becomes propor- 
tional to Jn+a + x (c7])l(cr]) a+ ^ . This follows from the integral representation (2) in which \p an (c,t) 
in the interval— 1 ^ t ^ 1 becomes proportional to ^(77), from the fact [1(20)] that 



/ 



1 V2^i'T(n + 2a+l)J ll+a+i (c v ) 
^(1-^(04- ^ £-— (8) 



for « > — 1. 

(d) The functions i//a/j(c,r/) are real for real 17, have exactly ft zeros in the interval (— 1, 1), 
and are even or odd functions of 77 according as n is even or odd. This is implied (not proved) by 
the limiting case c— >0 in which the limiting functions, the Gegenbauer functions of (4), are known 
from Sturm-Liouville theory to be real, to have exactly n zeros in (—1, 1), and to be even or odd 
functions of 77 according as n is even or odd. 

(e) The function i// a «(c, 77) cannot have a zero at ij=±l. This follows from the differential 
equation (1) and the entirety of ifj a n (c, 77) . For if it were to have a zero then, from (1), the first deriva- 
tive also must be zero there. And after I successive differentiations of (1) every /-h 1th derivative 
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must also be zero there, which would require that the function i// a «(c, 77) be identically zero every- 
where in 7). 

(f) The eigenvalues 6« w (c) of (1) are real, positive, and ordered such that b a o < b a \ < b a 2< • • ., 
provided that they are distinct. This is implied (again not proved) by the limiting case c— > 0, being 
given in the limit by (3). 

(g) The functions i/f« w (c, 77) behave asymptotically with real 77— > + oo as 



</faw(c, 17) = 



i w 2 tt+i r(a+l)^ ffll (c, 1) 



rc + a+1 

COS I CTj " 77 



V a n{c) 



(C7]) a ^ 



1 + 



m 



(9) 



This follows from the known asymptotic behavior [5] of integrals of the form (2). 

(h) The spheriodal wave functions S m , m+n(c, 7}) and the periodic Mathieu functions Se n (c, 7]) and 
So w + i(c, 7)) are special cases related to i//« w (c, 17) as follows 2 : 



^aw(c, 7)) = 



^m, m + n \C> J}) 
( 1 _ 7? 2)m/2 ' 

Se w (e, 7}), 

SOn+i(c, 7]) 

{ (i-,i)i/f 



a = m=0, 1, 2, 3, . 



«— - 



(10) 



1 

2* 



This is evident from the integral equation for the spheroidal wave functions [6(4)] and for the 
periodic Mathieu functions [6(11, 12)]. 

3. Construction of the Functions 

On the interval (— 1, 1) the Gegenbauer functions T%(r}) form a complete set with respect 
to all functions that are square-integrable with weight factor (1 — 7) 2 ) a , and are orthogonal with 
weight factor ( 1 — r? 2 )" [1(6)]: 



{ 



., (wmwr?(^= gtt 2 £ ) 1 *W 



(11) 



They are even or odd functions according as k is even or odd. Hence the spheroidal functions 
*/f«n(c, tj) can be expanded in Gegenbauer functions on ( — 1, 1) [1(8)], 



iffan(c,7])= ]£ d k (c\an)T%rj), 

k=0, 1 



(12) 



where the prime denotes summation over only even or odd values of k according as n is even or 
odd. The function i//«n(c, r\) becomes proportional to Tg(r)) as c— »0, hence 



2 The connection between the notation of Meixner and Schafke [3] and that used here is as follows: 



6 aM (c) = A£ + „(c 2 )-c*(a+l)+c 2 , <//„/,(r, 1?) a (1 --n 2 )"" l2 ps~^ n (r?; c 2 ) ; 

6 -i»( c ) = X »(-r) + T' i//-|,,,(c,7?)°c C e,,(cos -»ij; H; 



61 «(c)-A-»- 



1+ — . 
4/ 2 



^,»(c, i?) ^! — r} 2 ) 2 se„ + ,(cos x ^~t)' 
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d k (c\an) -^0,k^n. (13) 

Expansion (12) converges rapidly on the interval ( — 1, 1), but beyond that interval convergence 
becomes slower and slower with increasing 77. An alternative expansion that represents the func- 
tions uniformly on ( — °°, °°), by an argument used later for (40), can be obtained from (12) in the 
integral equation (2) by integrating term by term and making use of (8), 

i/j« n (c, 7]) = -— X a k (c\an) — , (14) 

V*n{C),£g fl (c7]) a+ 2 

where the expansion coefficients a k are related to the dk by 

, , x . fr T(A- + 2a+l) 7 , , x 
a k (c\an) = i k - n -+ - d k (c\an). (15) 



The coefficients a k are all real because the dk are real (the factor i k n is always real since k — n 
is always even). As c— >0, 

a k (c\an) -> 0, k ^ n, ( 16 ) 

by (13). For the special cases of a=m = 0, 1, 2, 3, . . . the coefficients a k - will agree with those 
tabulated by Stratton et al. [7] for the prolate spheroidal wave functions provided that the normaliza- 
tion is chosen to be such that 

J' i k - n a k (c\an) = l. (17) 

A-=0, 1 



With the right hand side of (17) replaced by V tt/2 the coefficients a a- for the special cases of 
a = ~i and +i are related to Blanch's [8] Mathieu function coefficients De^f and Do^liK 
respectively, by 



a k (c\-i,n) =i k - y|fl4 n) (c 2 ), 



a k (c\+h n) =»*-« y| (A-+1) DofcfHc*). 



(18) 



In addition to the increasingly poor convergence of (12) outside of ( — 1, 1) it requires special 
consideration for the case of a = — \ because T^ (r;) is infinite for k = 0. Hence the expansion 
(14) will be used almost exclusively hereafter. 

A recursion relationship [1(9)] for determining the a^'s in (14) and the eigenvalues b a n(c) 
is obtained from the differential equation (1), 



A k+2 a k+2 + (ban—Bk)ak+Ck-2Clk-2 = 0, (19) 
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where 

A-(A--l) 2 

A ' (2A + 2a-l)(2A + 2«+l) C ' 

« t^j-9 i n i 2A 2 + 2A(2g+l)+2g-l 

(A+2a+l)(A+2o: + 2) 2 
A ' (2A + 2a+l)(2A + 2a+3) ° ' 

It can be rewritten as a continued fraction for the ratio of two succeeding values of ak in two 
different ways, one in terms of increasing subscript, 

ttfc-2 , d i a fl fc+2 ' 

a* 
and the other in terms of decreasing subscript, 

Qfc-2 = ^ 

0«/i — Dk-2 T (>A--4 

Each must be the reciprocal of the other, which leads to an infinite set of eigenvalues for b ( , n . The 
best method now available for computing these eigenvalues, and simultaneously the expansion 
coefficients «a, is an iterative procedure that was first used by Bouwkamp [9(48)] for the case 
a = and later invented independently by Blanch |1()| for the Mathieu functions. All of the co- 
efficients ax in (19) must vanish for negative A, according to a result of Gegenbauer's [11, Sec 16.13 1 
in which he showed that any analytic function can be represented within its domain of analyticity 
by a series of the form (14). And for A =0 or 1 the //a , B*, C* in (19) become indeterminate at 
a = ±i. Therefore (19) should be supplemented by 

a k = for £<(), 

(22) 

r 2 3c 2 2^-h2 

#o— ^,^ = 2a + 2 + -^,C = f^c 2 , 
2a -to La+ 5 2a -f- 3 

the latter of which resolves the indeterminacy. For the special case of the Mathieu functions 

c 2 c 2 

(a = ±i) as computed by Chu and Stratton [2(100) J, where A k = C k = — and B h ■= A 2 -f— , the two 

computational procedures are quite different although the numerical results have been found to 

be exactly the same. The advantage of the new procedure here is that all orders a > — 1 are treated 

exactly alike. 

Once the eigenvalues are known the ratio — ^— for all k is also known, either from (20) or (21). 

dk-2 

The a a 's can then be computed from these ratios in terms of ao or a i for the even and odd spheroidal 

functions, respectively. The particular value chosen for ao or ai determines the scale factor by 

which the i//«/?(c, r\) functions are normalized. But no common agreement exists among authors 

on the choice of normalization factors. For this reason all formulas in this paper will be shown with 

arbitrary normalization factors that can be chosen at will. For those numerical results shown whose 

values depend upon a choice of scale factor the ifj (Y n(c, 7]) functions will be normalized to unity 

over the interval (— °°, °°). 

The eigenvalues v a n(c) can be computed by equating the two expressions (12) and (14), or 

their derivatives, at n = 0. As rj approaches zero the Bessel functions in (14) cause all terms to 

vanish except for k = Q. Hence for even rc = 2r, 
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, (-iy V^a (c\a, 2r) /oax 

""• 2r(C ' 2«r(a+i)*„,„(c,0)" (2l5) 

For odd rc = 2r+l the t//«»(c.T/) vanish at 17 = 0, so (14) must be differentiated with respect to 17 
before letting 17 approach zero. All terms vanish except for k=l, hence 

, v_ . (-l) r V^ai(c|q, 2r+l) 
" a , 2 r + ,(c)- lC 2~T(a+*)*;,, r+1 (c,0)' (24) 

The «^«n(c, 17) functions and their derivatives at 17 = are found from (12) to be 

*-• 2r(c ' 0) ~ wz 2 rTTT^TTT a2 ' (c|a ' 2r) ■ (25) 



2-V£«r(/+a+i) 



(26) 



Also, 



^«(c) = i» , |v«(c)|. (27) 

4. Double Orthogonality 

An important property that will be established here is that the spheroidal functions i//«/i(c, r/) 
for any order a > — 1 are orthogonal with weight factor |1 —r] 2 \ a on two different intervals (— L, 1) 
and (— °°, °°) simultaneously. This double orthogonality property was first recognized for the case 
of the prolate spheroidal wave functions of order zero by Slepian and Pollak [12]. A later investiga- 
tion by the author [6| revealed that this property was possessed by all of the prolate spheroidal 
wave functions of order m = 0, 1, 2, 3, . . . and by the periodic Mathieu functions. 

The differential equation for the functions S a , a +n(c, tj) defined by 

ba. g + njc, 7]) 
<ltan(c, V)=-Jl_ r] 2)al2 (28) 

\ 

is found from (1) to be of the singular Sturm-Liouville form 



[(i^ 2 )s;, a+ „(M)]'+ 



a 2 



b an + a(a+l) -c 2 T7 2 --T3— S«,<*+n(c, t?) =0. (29) 



Multiplication by S* a+p (c, 77) and subtraction of the product equation obtained when S a ,a+n(c, 17) 
and S a , e*+p(c, tj) are interchanged gives 

/ L V J ^7 / W a, a+i^ a, a+p *^«, a+H^> a. a+p) J V^«/> Oa/7 )»^ a, a+/?»J «, «+/>• (p") 

Integrating with respect to 17 over a real interval and changing from S«, «+«(c, 19) back to i//««(c, 17) , 

(l-77*)|l-77 2 KtK»(c, T?)<MC, 1?)-*«(C, TJ)^(C, I))]!?!,, 

C r) 2 
= (b a p — ban) tyan (c, 7j)l/fap (c, T/) 1 1 — 7} 2 \ a dr). (31) 

The functions ifj a n(c, rj) will be orthogonal over any real interval (171, r)>) for which the left hand 
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side of (31) is zero, provided that the b atl are distinct. Distinctness of the b a n can usually be deter- 
mined numerically. 

(a) Orthogonality on (— 1 , 1 ) 



When T7 1 = — 1 and 172= 1 the left hand side of (31) vanishes for all a> — ] because 
(1 — i7 2 )|l — 7j 2 | a vanishes at both end points and r/i a n(c, 17) and its derivative are finite. Hence, 
if the b all are distinct, 



/; 



*«(c f *?)<Mc, 'n)(l-7? 2 )«d'n = 0, n^p; (32) 



i.e., the i// a/ ,(c, r\) functions are orthogonal on (— 1, 1) with weight factor (1 — r) 2 ) a . 
(b) Orthogonality on (— °°, °°) 

From the asymptotic behavior (9) of i// m ,(c, 77) a simple calculation shows that 



(1 -T7 2 )|l— rf\ a Uf' nn {c, 7])\\i a „{c, 7))-\\) an {c, 7)) #ip(c, J])] 

\Van(c)\\v ap (c)\c^ +i 



sin — - — 7T + () — 

2 \7) 



(33) 



as 17— > + °° 9 7] real. When p — n is even the first term is zero for all 17, leaving 

(l-r ? 2 )|l-T ? 2 |-[^(c, 7))$ap(c, T7)-l// a „(c, l))^ ap (c, 7?)] = 0(1/7)). 

The left hand side of (31) then vanishes for 171 = 1 and 1)2 = °°, hence 

I i//«„(c, T?)i/fop(c, 17) 1 1 — 7} 2 | a dr) = 0, p — n even and 7^ 0. 



(34) 



Similarly on the interval (— °°, — 1), because of symmetry of the functions. Combining the three 
intervals (—00, —1), (— 1, 1), and (1, 0°), 



I <Jf«n(c, r))\\) al >{c, t?)|1 — r) 2 \ a din = 0, p 7^ n, 



(35) 



including p — n odd because i^« M (c, T7)i//«,,(c, 17) is then an odd function on (— °°, °°). Thus the func- 
tions i//«„(c, 17) are orthogonal with weight factor |l — r) 2 \ a on the interval (— °°, °°), provided only 
that the b att are distinct. 

5. Completeness 

(a) Completeness on (— 1, 1) 

The functions S a , a +n(c, t?) are complete on the interval (— 1, 1) with respect to all square- 
integrable functions. This is a property [13 1 of all Sturm-Liouville systems of the form 

(pn')'-h (kw — q)n = 0, puu'\ (t = puu'\t, = 0, (36) 
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and hence applies to the differential equation (29) and the boundary condition that S a ,a+n(c, n) 
must remain finite at the singular points f) =±1. Completeness of the S a ,«+n(c, 77) functions 
implies completeness of the i//a W (c, r\) functions. 

(b) Completeness on (— 00, 00) 

The fact that the functions i//aw(c, 77) are complete on (—1,1) with respect to the class of 
functions that are square-in tegrable with weight factor (1 — rf) a will now be used to show that 
they are also complete on (— 0°, 00) with respect to all functions of the form 



Fa(cri) = J e icr > ( (\ -t 2 ) a f(t)dt 



(37) 



for a > — I. The only restriction on f(t) is that it be square-integrable on (—1,1) with weight 
factor (1 — t 2 ) a . It can be expanded in a series that converges in the mean with weight factor 

(i-« 2 )°, 



/here 



/00 = y ^-^(c,i), -i<t<i, 

bn(c) = -^44 f /(*)*»(C, 0(1 - t 2 ) a dt. 



(38) 
(39) 



Hence for all real values of 17 the absolute value 



F a (cri) ~ J M»„n(c, tj) = f ' e^'(l - t 2 ) a f/(0 - 2 — <Mc 

„-0 J-\ L » =0 ^w 



/: 



(l-t 2 ) a 



dt 



bn 



n=0 ^<* n 



A (40) 



is exactly zero. It is concluded, therefore, that 

F a (C7)) = 2 b„(c)\p a „(c, Yj), — co < tj < cc, 



(41) 



where the series converges to F a (cr)) not just in the mean but uniformly on (— 0°, 0°). Thus the 
functions i//«n(c, 17) are complete on (— 0°, 00) with respect to all functions F a {cr\) of the form (37). 

6. Normalization Factors 

The double orthogonality property of the spheroidal functions results in two different sets of 
normalization factors, one on (—1,1) and another on (— 00, 0°). Both will be developed here. 

(a) Normalization factors on (— 1, 1). 

The orthonormality relation on (—1,1) will be designated by 



/: 



l/*cm(c, 7))t|»ap(c, 7)) (1 — rf) a dr) = A tt «(c)8„ 



(42) 



where A an (c) denotes the normalization factor for the spheroidal function of order a and index n. 
These normalization factors can be computed directly from the integral (42) for each p = n by 
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expanding i//«„(c, 77) in Gegenbauer functions as in (12) and making use of their orthogonality 
relation (11), obtaining 

A "" (C)= t (k+a + Dllk^a+l) ^^' (43) 

/»• = (), 1 V V ' 

where the a//s are the expansion coefficients (15). 

(b) Normalization factors on (— °°, °°). 

Each of the normalization factors on (— °°, °°) can be derived in terms of the corresponding 
normalization factor A ttW (c) on (— 1, 1) by making use of the double orthogonality property already 
established. The orthogonality integral (35) can be written 

i//«„(c, T))\\f ap (c, r?)|l -7] 2 \ a dr] 
J -00 

= I' \— Vt f " e ja *\l—n*\ a 1'<m(c 9 v)dv] (l~t 2 ) a $ap(c, t)dt, (44) 

J \ LVapKC) J -oc J 

by replacing \fj ap (c, rj) by its integral representation (2) and interchanging the order of integration. 

From the orthogonality property (35) on (— °°, °°) it is clear that this must be zero for all p 9^ n. 

And from the other orthogonality property (32) on (— 1, 1) it follows that the bracketed expression 

in (44) must be orthogonal to every i// a/ , (c, i) except i// a/ ,(c, t) . Then, from completeness on (— 1, 1) , 

it is concluded that the bracketed expression for p = n must be proportional to i//«„(c, t) everywhere 

on (-1, 1): 

1 f x 
— e ic *l-T) 2 hMc, T7)drj = y„ w (c)iMc, t), -1<*<1, (45) 

Van(C) J-x 

where y a n(c) is the proportionality constant.'* Thus the orthogonality relation on (— °°, °°) is 

tyanic, r))$ ap (c, v)\l-V 2 \ a dv = yan(c)Aan(c)Snp; (46) 

J — DO 

that is, each normalization factor on (— 00, °°) i s proportional to the corresponding normalization 
factor on (— 1, 1), the proportionality factor being the real positive number y a n(c) defined by (45). 
The new numbers y a n(c) representing the ratio of the two normalization factors are charac- 
teristic of the functions i// a »(c, 17), much like the eigenvalues b a n(c) and v a n(c), even though their 
defining equation (45) is not an integral equation in any ordinary sense. However they are by far 
the most difficult of the three to evaluate numerically for arbitrary a > — 1. In a few special cases, 
notably a = — i, 0, i, and 1, it becomes practical to evaluate them by evaluating the normalization 
integral (46). Such a method was actually used to obtain the first published values for the case of 
the prolate spheroidal wave functions of order zero [12, p. 59] (they were not recognized there as 
being a new kind of characteristic numbers because for a = they are not, being related simply 
to the eigenvalues of integral equation (2)) and the Mathieu functions [6 (60, 61)]. But for arbitrary 
a the only practical method to evaluate them appears to be from their definition (45), 



e"' T "|l -r) 2 \ a $an(c, 7))<h) 



yan(c) = J — 7-TTT~a ~ - 1 < / < 1 , (47) 



3 An early draft of the present paper was communicated privately to J. Meixner, who became sufficiently concerned about the lack of convergence of the integral 
in (45) for a 5 1 that he looked into the possibility of an alternative derivation. From some unpublished notes that he kindly sent me after this paper was submitted 
for publication it appears that he has succeeded in deriving the relationship (45) by an entirely different approach. He has given the integral from — °oto— 1 and 
from 1 to oo a meaning by assuming that t has a negative or positive imaginary part, respectively, and then approaching the limit of real t, thereby avoiding the difficult 
problem of proving summability. 
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at some convenient value of t in the interval (— 1, 1). In particular, at t = this reduces to 



■y«,2r(c) 



2 I" (7> 2 -1 )«.//„,. > r (c, 7,)<fy 

l=-i! — -— 



J'«,2r(c)^a,2r(c, 0) 



(48) 



for even ai = 2r. For odd n — 2r -f 1 the ratio (47) becomes indeterminate at £ = because 
i/>a,2> + i(c, 7)) is an odd function of r\. It can be made determinate by applying L'Hospital's rule 
at £ = 0, 



i2c T7(l7 2 — l) a l//a,2r+l(c, 7]) dr} 
Ya,2r+l(c) - 1 = ! — -^ J 777 



^a,2r+l(c)l//^ 2r+1 (c, 0) 



(49) 



for n = 2r + 1. The problem of evaluating y a /j(c) thus reduces to evaluation of the integrals in (48) 
and (49). 

The above method for determining y a n(c) was introduced originally by the author [6(42)] for 
the special case of the spheroidal wave functions of integral order a = m. It has been confirmed 
and generalized recently in an elegant and quite different way by Meixner [14], who showed 
that for integral a = m an expression equivalent to (45) holds for the most general solutions of 
the differential equation (1) for arbitrary complex values of c, of b, and of the characteristic exponent 
v in place of the integer n. 

The integral in (48) does not converge at infinity for a ^ 1, nor does the one in (49) for a ^ 0, 
but they both appear to be Cesaro summable (C, m + 1) for all — 1 < a ^ m because of the 
sinusoidal asymptotic behavior (9) of \p an (c, -n). They can be evaluated by expanding i// ttM (c, 77) in 
Bessel functions by (14) and integrating term by term, 



J*,2r{c) — 1 ■=- 



2 V27 



|<2r(c)l//a,2r(c, 0)|/-0 



2 a 2/ ,(c|a, 2r)h, 2} ,(c), 



(50) 



7a, 2r+l (c) - 1 
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K,2r+l(c)</>a,2r4-l(c, °) I ;^0 

where the integrals I a k(c) and K a k(c) are defined by 



Wc)=f(V-l)«%^ fl 



£ a, l> + l (c\a, 2r+ \)K a , 2l ,(c) 



K*(c) 



( CT? ) a + i 

It 2 -\\a J k + 2 + *-i (O?) , 

J, ^ ' } <«)-* ^ 



(51) 



(52) 
(53) 



For k ;= the integral / u o(c) in the range — 1 < a < can be evaluated by replacing the Bessel 
function in (52) by its representation [11, p. 48(2)] as an integral on (0, 1) and interchanging the 
order of integration, 



-a + i 



lav(c) r(a + i)r(i)J„ 



(i -t 2 ) a 



\: 



(tj 2 — l) a cos cyjtdiq 



dt. 



(54) 



The inner integral on (1, °°) i s an integral representation for the Neumann function [11, p. 170(4)], 
hence 
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'-(«)— fa-' 1 )- ^iSk 



(ct) 



(It 



cos avr 



S,nOT | (1 " f) (c*)« + * dt 

Jo (Ct)" + i 



(55) 



the latter expression coming from the definition of the Neumann function [11, p. 64(1)]. Expanding 
the Bessel functions in power series about the origin and integrating term by term gives 



i < x r(q+n ^ (-i)' +l 



2« + i '■ cos a7r . 



/! 



sin a;rr(/-a) /c\- 2 "-' !'(/ + ?) 



flr(/-a + t) \2 



P(Z+a + f) 



(56) 



The integrals I a k(c) and K a k(c) for all A: can now be obtained from I a o(c) by recursion in k. By 
use of the recursion relationships among Bessel functions [11, p. 45(1,2)] the following recursion 
relationships among / q a(c) and Kak(c) can be established, valid for all a > — 1: 



K a , k - i (c)+K ak (c) = (2k + 2a+l)U(c), 

Ka,*-2(C) -Kak(c) = (2a + 1) I ak (c) + 2c/^(c) , 

£«a(c) = klak{c) —cl'akic), 
{k+2a+3)l a , k+ 2(c) + cl' a , k+2 (c)=kl ak (c)-cl' ak (c). 



(57) 
(58) 
(59) 
(60) 



The last relation (60) is a differential equation for I a ,k + 2(c) in terms of I a k(c). Its solution is 

/a, fc+ 2(c)= Zg c , + ^ + 3 j o t* +fc+ «/ alk (0*-/*(c) i (61) 

Repeated application of (61) to I a o(c) in (56) suggests the appropriate form of the power series 
expansion for I a ,2p(c) , which is then easily shown by induction to hold for all nonnegative integers p. 
Itsyth derivative, obtained by differentiating term by term with respect to c, is 

_ r(q+l) ^ (-iy^ \ S ma7rr(l-a)(-l + a + i) p (2l-2a-j) j /c\- i "-' 
i a ,2 P yc> 2 >+«+» cos «r £ /! L (/ + p)!r(/-a + i) 



r(/+l)(-0,(2/+i-i) j 

r(/ + a + |)r(/ + a + |+p) 



c\ 2 '-J 

2, 



(62) 



where (z)p = z(jz: + l) (z + 2) . . . (z~\-p — 1), (z)o— 1, is Pochhammer's symbol. I a ,2p(c) and 
K a ,2p(c) for — 1 < a <0 can thus be computed by (62) and (59), respectively. 

To obtain I a ,2p(c) and K a ,2 P (c) for orders beyond the region of integrability — 1 < a < of 
(55) the following recursion relations are valid for all a > 0, 



/«*(C) = 



1 



c 2 2/c + 2a + l 



/a-l,fc(c) ~ 2A , + 2a+ j- /a-l,fc+2(c) — /i_ lfjb (c), (63) 



JC*(c)«*£^/._^ 



(64) 



the first being obtained from the definition (52) and the Bessel function recursion relations and the 
second from (59). It has been found by direct substitution that the series (62) satisfies the recursion 
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388-513 O - 70 - 6 



relation (63). From this it is concluded that (62) is valid for all a> — \, even though the integral (55) 
from which it was derived contains a nonintegrable singularity at t = for all nonintegral values 
of a greater than zero. As an indication of the rapidity of convergence of (62) and of (50) and (51), 
and hence of the practicality of the formulas developed here, the series (62) always was found to 
converge to 16 significant figures in less than 23 terms for the computations shown later, and the 
series (50) and (51) to converge to 16 figures in less than 19 terms. 

The series (62) becomes indeterminate when a is an integer or a half-integer. The indetermi- 
nancy at any of the integers a = m can be resolved by taking the limit as a approaches m. For 
ct = 0, for example, the limiting value /oo(c) is 



1 



/oo(c) 
in agreement with the power series expansion of 



c 21 



\2 Vtt£ V ; (2/+l)(2/+l)!' 



(65) 



1 2 (it q . 

C V '77 VZ 



obtained by direct evaluation of the integral (52). At a = — \ the indeterminancy again can be 
resolved by taking the limit, obtaining 



l -^-ir,%y vl W\ m+l 



)-*(/ + « -21n- 



2/ 



(66) 



in terms of the Digamma function i//(z), which is identical to the power series expansion [11, 
p. 150 (1)] of 

~2 J«(2) Yo \2 

obtained by direct evaluation of the integral (52). But in the special cases of integral and half- 
integral orders, and these cases only, there are more direct means available for evaluating y a n(c) 
than by the series (50) and (51). For all integral orders the numbers y a n(c) are found to be express- 
ible quite simply in terms of the spheroidal eigenvalues b a n(c) of the differential equation and 
Van(c) of the integral equation [6(23)]. For the half-integral orders a = ±i the numbers y a n(c) 
have been found [6(31)] in terms of Blanch's joining factors f e , w (c 2 ) and f , n +i (c 2 ) to be simply 



yan(c) 



l+f e , n (c 2 ), a = -i 

l+fo,n + l(c 2 ), a = + i 



(67) 



From Blanch's expansion of/ e , n (c 2 ) and f ,n+i(c 2 ) in terms of Bessel functions it can be shown 
that (67) is identical to that obtained from (50) and (51). The integrals I a k{c) and K a k(c) can be 
expressed in closed form for half-integral values of a, starting from I-i,2p(c) [15(6.552.6)] 



/-i2p(c)=-|j P (|)y P (0, 

from which /£_! ,2p(c) is obtained by recursion relation (59), 



K-. 



(0=-f[; p (|)y p+I (|) + y p+1 (|)F p (|) 



(68) 



(69) 
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and from which /*, 2p(c) and K i zp(c) are obtained by recursion relations (63) and (59), 
respectively, 



h.2p(c) 



4(2/;+ I) L "\V "V2 



4^(9 + ./,,©>,,({ 



2p±l 
c 



„(|)^,(|) + y„,(f)v,.(|))]. ao, 



K|.»(c)-f 



./,.©v„(5) + 7„.,(0t„.^ 



wiK'iD+m) 1 ^))} ™ 



7. An Extremal Property 

An important property of the prolate spheroidal wave functions of order zero was established 
by Landau and Pollak [16, Th. 1], namely, that the first /V-f 1 of them are the /V+ 1 linearly inde- 
pendent functions of the form Fo(cr)) that are most concentrated in the interval (— 1, 1 ). By the 
term "concentration" of a function of the form F a (cr)) defined by (37) will be meant the value 
of the ratio R a (c) , 

[' |JF«(cii)|"(l-i|«)^i) 

Ra{c)=j=* (72) 

J |F (t (c7,)| 2 |l-T, 2 hiT, 

It will now be shown, by quite a different and much more simple argument, that for any order 
a > — 1 the first /V+l of the spheroidal functions i/f« w (c, tj) are the N + 1 linearly independent 
functions of the form F a (cr\) that are most concentrated in the interval (— 1, 1). 

After expanding F a (cy]) in spheroidal functions as in (41) and making use of the two orthog- 
onality relations (42) and (46) the concentration ratio (72) becomes 

PC 

2 b n b*Aan 
Ra = ^ — (73) 

2 bnbtyanAan 
n = Q 

It is determined entirely by the expansion coefficients b n of the function F a {crj). Hence a neces- 
sary condition for it to attain its largest value is that it be an extremum with respect to the real 
and the imaginary parts of each and every one of the bn$; i.e., since b n and bn are linearly related 
to the real and imaginary parts of b n , 

dR a dR a n A 1 o O /r^x 

^ = ^* = ' P = 0,l,2,3, ... . (74) 

Rewriting (73) in the form 

|wi(ft a r»»-l)A„„ = 0, (75) 

(1=0 

then differentiating through with respect to b,, or b * and imposing condition (74), 

b* ( Raja,, - 1 ) Aap = b P {R a y ap - 1 ) A„„ = 0, p = 0, 1 , 2 , 3, . . . . (76) 
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The normalization factors A a/ ,(c) are never zero, hence either 

b p = b*=0, (77) 

or 

Rajap- 1 = 0. (78) 

The second condition (78) can be true for at most only one value of p, in which case condition (74) 
for all other values of p requires that (77) be true. It is clear, then, that an extremum of R a (c) can 
occur when, and only when, F a (cr)) is proportional to one of the spheroidal functions i//«n(c, 77). 
The value of R a (c) at each of these extrema is, from (78), 

K«ext(c)=7nvv " = 0>1*2,3, . . .. (79) 

Jan \C) 

The largest possible value of R a (c) is that for which y a n(c) is least, which is y a o(c). Hence the 
spheroidal function for n = is the single most concentrated function of the form F a (ciq). 

By the term "linear independence" of a set of functions is meant that they are mutually orthog- 
onal. Thus to find the second most concentrated function of the form F a (cr)) that is linearly 
independent of i//«o(c, 77) on (— °°, °°) means to find the most concentrated function of the form 
F a {ciq) that is orthogonal to i//«o(c, 77) on (— °°, °°). From the orthogonality of the functions i// a »(c, 77) 
on (— °°, °°) and from their completeness with respect to functions of the form F a (cr)) it follows 
that all functions of the form F a (cr)) that are linearly independent of i// a o(c, 77) can be expanded in 
the set of i// tt „ (c, 77) functions that remain after i//« (c, 77) has been removed; i.e., after setting bo = 
in (41). Repeating the maximization process above leads again to conditions (77) and (78) but 
restricted now to p=l, 2, 3 . . . . The largest possible value for R a (c) is again that for which 
y an (c) is least, which now is y«i(c). Hence the spheroidal functions for n = and 1 are the most 
concentrated pair of linearly independent functions of the form F a {cr\). By this process of sys- 
tematically depleting the set of functions i//« w (c, 77) it is concluded that the first N+ 1 of them are 
indeed the 7V+ 1 linearly independent functions of the form Y a{cr\) that are most concentrated 
in the interval ( — 1, 1). 

8. Spheroidal Identities 

All of the identities to be derived here will be obtained from the integral equation (2) and the 
following property of completeness of the functions i//cm(c, 77) on (— 1, 1): 

8(l - o = (i - ,*)<« f fe&iilfakil (80) 

for — 1 < t < 1 and — 1 < t' < 1, obtained by expanding the delta function 8(t — t') in the ortho- 
normal functions *//««( c, £')/VA«„(c) with expansion coefficients (1 — t 2 ) a \\} an {c, t)ly/A an (c). 
Although such an expansion is open to question on grounds that the delta function is not square- 
integrable it does appear to be valid as the limiting case of some square-integrable function that 
approaches the delta function. 

The principal identity from which the others will be derived is obtained by multiplying the 
integral equation (2) for argument 77 by its complex conjugate for some other argument y, then 
dividing by the normalization factor A an (c) , and finally summing over all n, 



f ¥^# fc»(c, 7,)^»(c, y) 

„ = tlanyC) 



/:/: 
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(l— t 2)a y *Pan(c, t)\\) an {c, t' ) 



Aan(c) 



dtdt'. (81) 



The bracketed expression is 8(t - /'), by (80), which reduces the right hand side to an integral 
representation of a Bessel function [11, p. 48(4) |. Hence one obtains the identity 

S W" '""" <r - T ' H '"" (r ' V) ( ° U> (e[T,-r])^ (82) 

for — oo<^<oo i — 00<y<00. 

An important special case of (82) occurs when y=rj, giving 

„4j a«,(c) "^'^ r(«+i) ( 83 > 

for — =° < 17 < oo. At 17 = the product v a n(c)tjfan(c, 0) for odd n is zero and for even n is propor- 
tional to «o, by (23), hence 

» al(c\a, 2r) 

X~A 7ZT = i r ( Za + 2 )- (84) 

This identity provides a means for an overall numerical check on computed values of the normaliza- 
tion factors Aon(c) alone for even n. To obtain the corresponding identity for odd n construct a 
new identity like (83), 

* ka„(c)| 8 „ . r(a+i)r(i) 

</»aUC, 7j) = < ... (85) 



& a«(c) """ v -" i'(«+D 

for — rc<i7<=°. At 17 = this reduces, by (24), to 

» a?(c|a, 2r+ 1) 

2 A 7- r i =(a + i)r(2a+2), (86) 

which provides a means for checking numerical values of the A aw (c)'s alone for odd n. 

An identity for the eigenvalues v a n(c) alone can be obtained by multiplying (83) by (1 — T7 2 )" 
and integrating over ( — 1 , 1 ) , 

?„ Mc)l I r(«+i) j' (87) 



An identity for checking the new numbers y a n(c) for even n can be obtained from (82) by 
letting y=0, multiplying by (iq 2 — l) a , and integrating over (1, 00), 

- ( _ 1)r K2r(c)|>» 2 (c,0) [ya2r(c) . 1]=2 a + }r(a+l)rft)/„o(c), (88) 

^J A tt ,2r\CJ 

where / a o(c) is as defined by (52). For odd rc differentiate (82) before letting y=0 and then multiply 
by n(77 2 — 1 ) a and integrate over (1, °°), 

y ( -i)r + i l"-^'' c ^;'^' (c - 0) [ y „., 4l ( c )-i]=2-ir(«+i)r(l)r/;(0. (89) 

£51 A„,2/ + i(c) 
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Another identity can be constructed by multiplying integral equation (2) by i//«n(c, y)/V A an (c) 
and summing, 

2 ?"#*««(c, i})*«(c, y)=«fcw, (90) 

for — °°<T)< 00 , — 1 < y < 1. For y=0 this becomes, from (23) and (24), 

A a (c\a,2r) 2 T(a + 3) 

S (_1) A Q , 2r (c) ^-^c, V )— - m — , (91) 



r=0 



\ (-l) r \ /TT— *a f 2r+l(c,T7)- fTTx ^ 



(92) 



for — oo < 7] < oo. 

In the limiting case where c —> and 17 is such that 07 = £, 



^,(c^)^ r(/t + ^ +I) a^(r ? ), (93) 



( n -h a + f) 1 ( /i H- za + 1 ) 

1 . / Y \ 
\v an (c)\\lfan(c, T)) -> V^T fl„ " + "^ (95) 

Then (82) for c^ with ctj = £ and cy = z reduces to 

- (n+q + i)r(Ai + 2a+l) J,^U) 7„ + « + *U) = 2 g+ *r(q+l)rq) Ja+iU-z) 

which is one of Gegenbauer's addition theorems [11, p. 363(2)], and (90) for c— * with ct7 = £ 
reduces to 



277 J i»(n+a + i) • / " +n ; + y ) r-(y) = e*0r, (97) 

»=0 b 



which is another [11, p. 368(2)]. 



9. A Special Case 

For the special case in which c=— , a=l, and n=q — 1, with g=l, 2, 3, . . ., an exact 

closed-form representation can be obtained for the spheroidal function and each of its associated 
numbers. Its importance lies in the fact that it provides an independent means for an exact numeri- 
cal check. 

The spheroidal function sought can be obtained from the differential equation (1) and the 
condition of finiteness at 17 = ±1 by letting 

*Mc, 17) = 2 , (98) 
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which leads to 

i — J)- / 

When a= 1 this reduces to 

(l-^ 2 )tt? tt +(6in-+-2-c 2 i7 2 )ain = 0. 

</77 



For c = -^— and assuming for the moment that 



(100) becomes, for n = q — 1 , 



*...- ?-?» - 2 - 



(y) ki,«-i 



whose solution is 



M-l, r/-l 



hence 



«pi.g-i 



£77 

2 ' 



077 \ , COS 077 



COS #77 

sin T 



)dd 



1-TT 



(99) 



(100) 



(101) 



(102) 



(103) 



(104) 



where k q is a normalizing constant. Expression (104) satisfies the differential equation (1) and is 
finite at 7} = ± 1, so it is clearly a spheroidal function. Furthermore it has exactly q—\ zeros in the 
interval (— 1, 1), so the index of the spheroidal function must be q~\. Thus the validity of (101) 
is established. The spheroidal function (104) and its eigenvalue (101) were obtained earlier by 
Flammer [17]. 

The spheroidal function (104) satisfies also the integral equation (2), which can be verified 
by direct calculation of the integral, 





!> 


if-nt {1 _ t2] 


COS077 

sin 2 . 

i-i- *" 


qir 


cos 
sin 
1- 


qrr 
-r, 2 










(105) 


This is the integral 


equation in question, 


from which 


it is evident tl 


lat its 


eigenvalue is 












V\ ti 


- (?)=<■ 


i 4 

7-1 

#77 














(106) 


The normalizat 


ion factor 


on (-1, 1) 


in terms of 


the consta 


nt k q 


is found 1 


y direct 


cal 


leu 


lation 


of the integral (42) 


to be 


Ai,«- 


/<777\ k$ 

At)- 2 


Cin 2^77. 














(107) 
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And by direct calculation of the integral (46) one finds 

yi..-i (y)= 2 - < 108 ) 

The above results provide an interesting indication of the way in which the integral in the 
definition (45) of y an (c) can behave outside of the open interval — 1 < t < 1 . For q odd, for example, 
the left-hand side of (45) can be evaluated in closed form, 



(-i 



qir 



air 

COS „ * 



2 ' «=* 



= 2-^^ (-1) * |[8(*-l)+8(*+l)], (109) 

which agrees exactly with the right-hand side of (45) within the open interval — 1 < t < 1. It does 
not agree at the endpoints £ = ±1, however, as evidenced by the delta function singularities that 
appear there. 



10. Numerical Results 

The formulas developed above have been used to compute and check the spheroidal functions 
and their associated numbers for values of c up to 577. Representative numerical results for c = 6 
are shown below, all of which were computed on the Triangle Universities IBM 360/75 computer 
by Guy A. Myers at the North Carolina State University. 

The first nine of the spheroidal functions i// a/ ,(c, tj) are shown in figure 1 for a = — i, 0,i, 1,2, 
and 3. They were normalized on the interval (— -°°, °°) by choosing A an (c) to be lly an (c). The 
corresponding functions for a = with c = 0.5, 1, 2, and 4 were shown earlier by Slepian and 
Pollak [12]. The characteristic numbers fe aw (c), j> aM (c), and y an {c) as a function of a are shown 
in figures 2, 3, and 4, respectively. 

Numerical results for a few representative cases are shown in tables 1 and 2. Table 1 shows 
computed values of b an {c) , \v an {c) |, and i// tt /, (c, 0) for even n and \\}' an (c, 0) for odd n, for a=— £, 
0, and +| at c = 6, n = to 13. Table 2 shows computed values for the reciprocal of y a n(c) for 
a = — h~ 2,~~ i, 0, i,i, f, 1, and f, at c = 6, ra = to 13. The reciprocal of y an (c), rather than y a n(c) 
itself, is shown in order to illustrate the relative concentration (79) of the functions. 
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^n< c .n> 



fenfc.n) 



tntcn) 




tn (c .n) 



*-.« 



t«w 





lihim 

'I/' 



Figure 1. The spheroidal functions fora = — 1/2, 0, 1/2, 1, 2, 3, c = 6, n = to 8, M)j</< A„„(c)= l/y„„ (c). 
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Figure 2. The eigenvalues b an (c) of the differential equation (1) for — 0.92 =S a «S 1, c = 6. 



Hw c) | 




Figure 3. 77ie eigenvalues v an {c) of the integral equation (2) for —0.92 *S a ^ 1, c = 6. 



log(r an (0-l) 




-2 



-3 - 
-4 



-n=0 



Figure 4. 7^e new numbers y an (c) defined by (45) /or — 0.92 ^ a ^ 1, c = 6. 



206 



Table 1 



n 


ban(c) 


|*Uc)| 


«Mc o) 

W„(c,0) 






a = -\/2, c = 6 







0.5737585781+01 


0.1074447786 + 01 


0.1155255760 + 01 


1 


.1661329843 + 02 


.1212324184 + 01 


.3686243582 + 01 


2 


.2598284316 + 02 


.1469587054 + 01 


-6927129388 + 00 


3 


.3290367966 + 02 


.1770806514 + 01 


-.2666895609 + 01 


4 


.3816092638 + 02 


.1311549557 + 01 


.3176164384 + 00 


5 


.4509186607 + 02 


.4656092304 + 00 


.1273298762 + 01 


6 


.5521934384 + 02 


.1117504031+00 


-.6410649522-01 


7 


.6785668990 + 02 


.2310300898-01 


-.1308328538 + 00 


8 


.8264753135 + 02 


.4243576673-02 


.3284620306-02 


9 


.9950842599 + 02 


.6981239327-03 


.5920679880-02 


10 


.1184102166 + 03 


.1037825612-03 


-.9512187852-04 


11 


.1393381242 + 03 


.1405973769-04 


-.1623325826-03 


12 


.1622835823 + 03 


.1748761994-05 


.1808459022-05 


13 


.1872412953 + 03 


.2010009631-06 


.2990162450-05 






a = 0, c = 6 







0.5208269159 + 01 


0.1023276503 + 01 


0.1153430212 + 01 


1 


.1600044275 + 02 


.1021309607 + 01 


.3662488818 + 01 


2 


.2535647864 + 02 


.9922435547 + 00 


-.7045705538 + 00 


3 


.3320419949 + 02 


.8229938914 + 00 


-.2977793219 + 01 


4 


.4072019426 + 02 


.4659781026 + 00 


.2864309440 + 00 


5 


.4977371212 + 02 


.1693510360 + 00 


.7958373899 + 00 


6 


.6118075689 + 02 


.4524678973-01 


-.3152431639-01 


7 


.7485286652 + 02 


.9966212672-02 


-.6222320400-01 


8 


.9065115936 + 02 


.1897100693-02 


.1388242306-02 


9 


.1085154453 + 03 


.3192404741-03 


.2472548565-02 


10 


.1284188799 + 03 


.4820488557-04 


-.3606307265-04 


11 


.1503474434 + 03 


.6605196410-05 


-.6113574747-04 


12 


.1742930029 + 03 


.8286721621-06 


.6275570287-06 


13 


.2002505133 + 03 


.9588946460-07 


.1033409530-05 






a = + 3/4, c = 6 







0.4524868470 + 01 


0.9514693801+00 


0.1157683567 + 01 


1 


.1554416968 + 02 


.7960998515 + 00 


.3741186029 + 01 


2 


.2556039130 + 02 


.6000021716 + 00 


-.7319049126 + 00 


3 


.3502025958 + 02 


.3687113331+00 


-.2915597899 + 01 


4 


.4493678220 + 02 


.1690739803 + 00 


.1800330864 + 00 


5 


.5631678092 + 02 


.5762489345-01 


.3736554493 + 00 


6 


.6963816767 + 02 


.1546055800-01 


-.1125847285-01 


7 


.8500027646 + 02 


.3450023959-02 


-.2071854851-01 


8 


.1024006606 + 03 


.6633536963-03 


.3879838010-03 


9 


.1218275965 + 03 


.1124131264-03 


.6700408720-03 


10 


.1432726639 + 03 


.1705749626-04 


-.8491880927-05 


11 


.1667303958 + 03 


.2345366168-05 


-1414240391-04 


12 


.1921972059 + 03 


.2949700757-06 


.1289490027-06 


13 


.2196706817 + 03 


.3419275479-07 


.2099163320-06 
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Table 2 



Reciprocal of y a n(c) for c = 6 



n 


a = -3/4 


a = -l/2 


a = -l/4 





0.9998903116 + 00 


0.9998968961 + 00 


0.9999006633 + 00 


1 


.9953591572 + 00 


.9957735784 + 00 


.9959984111 + 00 


2 


.9251077487 + 00 


.9361280337 + 00 


.9403541758 + 00 


3 


.6617807432 + 00 


.6867880908 + 00 


.6773573623 + 00 


4 


.4329978267 + 00 


.3559391049 + 00 


.2783892137 + 00 


5 


.1631458857 + 00 


.8787799051-01 


.4890053549-01 


6 


.1812019930-01 


.8375304321-02 


.4009596286-02 


7 


.1085823478-02 


.4742740153-03 


.2109141235-03 


8 


.4704479809-04 


.1948721509-04 


.8155744837-05 


9 


.1566525909-05 


.6167222470-06 


.2444286528-06 


10 


.4138876101-07 


.1554061438-07 


.5862729444-08 


11 


.8891403564-09 


.3195449002-09 


.1152410426-09 


12 


.1584025115-10 


.5466463632-11 


.1891529381-11 


13 


.2378498905-12 


.7904581458-13 


.2632527573-13 


n 


a = 


a = +l/4 


a = +l/2 





0.9999018826 + 00 


0.9999006439 + 00 


0.9998968665 + 00 


1 


.9960616432 + 00 


.9959753979 + 00 


.9957367356 + 00 


2 


.9401733901 + 00 


.9363969465 + 00 


.9290353072 + 00 


3 


.6467919491+00 


.5993323516 + 00 


.5374082010 + 00 


4 


.2073492168 + 00 


.1474458621+00 


.1005292486 + 00 


5 


.2738716624-01 


.1530259540-01 


.8498117614-02 


6 


.1955000733-02 


.9617772918-03 


.4747110329-03 


7 


.9484876556-04 


.4292801980-04 


.1948802000-04 


8 


.3436783286-05 


.1454316917-05 


.6167231119-06 


9 


.9732115988-07 


.3886092897-07 


.1554061496-07 


10 


.2218980545-08 


.8416064934-09 


.3195449004-09 


11 


.4166226283-10 


.1508525403-10 


.5466463632-11 


12 


.6557478591-12 


.2276034041 - 12 


.7904581458-13 


13 


.8780377122-14 


.2931302844-14 


.9790506399-15 


n 


a = + 3/4 


a = +l 


a=+5/4 





0.9998902878 + 00 


0.9998804332 + 00 


0.9998665602 + 00 


1 


.9953272353 + 00 


.9947106203 + 00 


.9938277694 + 00 


2 


.9175535778 + 00 


.9009322432 + 00 


.8776555430 + 00 


3 


.4640663909 + 00 


.3839718982 + 00 


.3032540474 + 00 


4 


.6605749034-01 


.4205022091-01 


.2605099768-01 


5 


.4681647613-02 


.2555876781-02 


.1381915965-02 


6 


.2342324766-03 


.1152638318-03 


.5647610537-04 


7 


.8852297954-05 


.4016362284-05 


.1817713005-05 


8 


.2616675768-06 


.1109412441-06 


.4695584821-07 


9 


.6217024284-08 


.2485740919-08 


.9925660313-09 


10 


.1213549245-09 


.4606583810-10 


.1746778699-10 


11 


.1981169523-11 


.7177204946-12 


.2597757239-12 


12 


.2745434078-13 


.9531862327-14 


.3306789087-14 


13 


.3270100000-15 


.1091853480-15 


.3643090846-16 
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